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We report an investigation to establish the physical mechanisms responsible for decoherence in 
the generation of photon pairs from atomic ensembles, via the protocol of Duan et. al. for long 
distance quantum communication [Nature (London) 414, 413 (2001)] and present the experimental 
techniques necessary to properly control the process. We develop a theory to model in detail the 
decoherence process in experiments with magneto-optical traps. The inhomogeneous broadening of 
the ground state by the trap magnetic field is identified as the principal mechanism for decoherence. 
The theory includes the Zeeman structure of the atomic hyperfine levels used in the experiment, and 
the polarization of both excitation fields and detected photons. In conjunction with our theoretical 
analysis, we report a series of measurements to characterize and control the coherence time in our 
experimental setup. We use copropagating stimulated Raman spectroscopy to access directly the 
ground state energy distribution of the ensemble. These spectroscopic measurements allow us to 
switch off the trap magnetic field in a controlled way, optimizing the repetition rate for single-photon 
measurements. With the magnetic field off, we then measure nonclassical correlations for pairs of 
photons generated by the ensemble as a function of the storage time of the single collective atomic 
excitation. We report coherence times longer than 10 (is, corresponding to an increase of two orders 
of magnitude compared to previous results in cold ensembles. The coherence time is now two orders 
of magnitude longer than the duration of the excitation pulses. The comparison between these 
experimental results and the theory shows good agreement. Finally, we employ our theory to devise 
ways to improve the experiment by optical pumping to specific initial states. 



I. INTRODUCTION 

Quantum memory is a key resource for many quantum- 
information protocols. Usually it is associated with the 
basic requirements for quantum computation 0, Q , but 
in recent years also quantum communication protocols 
started to rely on it. The requirement of memory was in- 
troduced in quantum communication as part of the idea 
for quantum repeaters 0, 3 , a possible solution for the 
problem of quantum communication over long distances. 
In this case, memory is essential to increase the prob- 
ability of success of the chain of conditional steps that 
underlies the protocol, and makes feasible scalable quan- 
tum networks. 

A significant step toward the realization of the quan- 
tum repeater idea was a proposal by Duan, Lukin, Cirac, 
and Zoller (DLCZ) for its implementation using linear 
optics and atomic ensembles |5j- The DLCZ protocol is 
based on the generation of single photons by spontaneous 
Raman scattering in atomic ensembles |6| . The detection 
of a single photon in the forward propagating mode her- 
alds the presence of a single collective atomic excitation 
in the sample, due to a collective enhancement effect. 
This excitation can be stored for a time up to the co- 
herence time of the ground states of the atoms and then 
converted back into a light field. Entanglement of distant 
ensembles in the excitation number basis is generated by 
interference 0, and extended to longer distances by en- 
tanglement swapping The final pairs of ensembles, 
far apart, can then be used for entanglement-based quan- 
tum cryptography 0, , probabilistic quantum telepor- 
tation and violation of Bell inequality. This proposal has 
received much attention in the past two years and several 
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In this article, we analyze the decoherence processes 
present in the DLCZ protocol, and describe experiments 
to mitigate the problem. We construct a theory for the 
decoherence process in the photon-pair generation. Par- 
ticularly, our analysis concentrates in its implementation 
with cold atomic ensembles, but many results should 
also apply to studies with room-temperature ensembles 
in vapor cells. We propose various strategies to increase 
the system's coherence time, and introduce experimental 
techniques necessary for its characterization and control. 
We also report the first experimental steps in this direc- 
tion, with an increase of more than two order of magni- 
tude in the coherence time with respect to the previously 
reported works with cold atoms [Til llH ll H I m ITsI] . 

The coherence times reported up to now by the sev- 
eral groups working on the implementation of the pro- 
tocol are all shorter or of the order of a couple of mi- 
croseconds. Furthermore, for all experiments to date, 
the reported coherence times are of the order of the ex- 
citation pulses duration. However, for using this system 
as a quantum memory, it is important to obtain storage 
time much longer than the excitation pulses. Moreover, 
for the DLCZ protocol to become a viable alternative for 
long distance quantum communication, long coherence 
time is crucial and major efforts are required to increase 
it. The main goal of the present article is then to provide 
the initial steps in this direction, and to establish several 
techniques and ideas for the next steps. 

Only two types of systems have been em ploy ed in the 
experiments up to now: vapor cell s Il2t 1 13. 1 151 and 
cold atoms in magneto-optical traps |lll Il4 lla. Il7l Il8| . 
In both systems, however, the experiments have not 
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achieved yet their respective state-of-the-art coherence 
times. The vapor-cell studies, for example, did not em- 
ploy paraffin coated cells 0, |2£j]; the coherence times 
were effectively limited to the time the atoms take to dif- 
fuse out of the excitation region, which is of the order 
of microseconds. Recently, high fidelity atomic quantum 
memory of the state of a lig ht pulse was achieved with 
such paraffin coated cells j2l| with memory times of up to 
4 ms. Coherence times of tens of milliseconds, however, 
are commonly achieved in this system |2^| , and there are 
reports of coherence times as high as one second [l9j. 
The difference in these values is largely due to measure- 
ments of decay of different coherent processes 0. How 
the coherence required for the generation of photon pairs 
from atomic vapors will decay as the atoms collide with 
the walls of paraffin coated cells is still to be determined. 

The use of atomic traps to generate photon pairs for 
the DLCZ protocol has the advantage of providing a high 
density of atoms distributed in a small spectral region, 
due to the suppression of Doppler broadening by the cool- 
ing process. This allows the use of excitation laser pulses 
tuned closer to resonance, which requires much less power 
and makes it easier to filter the excitation pulses from the 
Ramam-scattered photons. However, atomic traps also 
introduce a different set of complications. In the case 
of the magneto-optical traps (MOT) used up to now, 
the magnetic field of the trap induces decoherence on 
a timescale of the order or smaller than a few hundreds 
nanoseconds Il7l Il8| . The first results with the MOT 
magnetic field off are reported in the present article, with 
coherence times on the order of 10 /is. As will be dis- 
cussed below in detail, a better nulling of the magnetic 
field combined with optical pumping to specific Zeeman 
levels might increase the coherence time, in a straightfor- 
ward way, to hundreds of microseconds. 

Further improvements with MOTs would face the 
problem of diffusion of atoms from the excitation region 
and, most troublesome, from the MOT itself. This prob- 
lem can in principle be mitigated by improved cooling 
techniques. However, along these lines, it would be dif- 
ficult to increase the coherence time above a couple of 
milliseconds. A possible solution then is to use an op- 
tical dipole trap to hold the atoms during the write- 
and-read process. Hyperfine coherence times of hun- 
dreds of milliseconds have already been observed in such 
traps Hill. 

In the following, Sees. ITT1 and IIVI are devoted to theo- 
retical results and Sec. Illll to associated experiments. In 
Sec. Ill Al we give a general introduction to the photon- 
pair generation process behind the DLCZ protocol. In 
Sec. Ill Bl we derive a theory for the probability of 
joint detection of these photons pairs generated from an 
atomic ensemble in a magneto-optical trap. This theory 
is a direct extension of a previous theoretical treatment 
reported in Ref.0, to which we added explicitly the read- 
ing process and the Zeeman structure of the levels. In 
this way, we are able to model the action of the magnetic 
field over the atoms, and to study the dependence of the 



correlations with the light polarization. 

Section IIIII describes an experimental investigation 
leading to the nulling of the magnetic field in the photon- 
pair correlation measurements, with the subsequent in- 
crease in the system coherence time and degree of cor- 
relation. In Sec. IIII Al we describe a series of Raman- 
spectroscopy experiments to characterize the system and 
optimize the process of zeroing the magnetic field. We de- 
termine the set of experimental conditions that result in 
a good compromise between atomic density and magnetic 
field cancellation, which we used in the correlation mea- 
surements. Section IlII Bl describes then measurements of 
nonclassical correlations for the photon pairs generated 
by the MOT. We compare results with magnetic field on 
and with magnetic field off. The magnetic field off mea- 
surements present a higher degree of correlation, and a 
hundred times larger coherence time. We compare the 
shape of the experimental curves with magnetic field on 
and off to our theory, obtaining good agreement. We also 
show how the two-photon wavepacket that describes the 
detailed temporal structure of the photon pair generation 
is modified by the magnetic field. 

Finally, based on the procedure for comparison be- 
tween theory and experiment described in Sec. IIII Bl we 
formulate in Sec. IIVI a proposal to improve our experi- 
mental signal. We suggest using a combination of optical 
pumping to a specific initial state and polarization of the 
light fields to increase both our detection efficiency and 
coherence time. Section [V] is dedicated to our conclu- 
sions. 



II. THEORY 

The basic theory for the DLCZ protocol is described 
in Refs. [f| and_p. The general idea of the protocol is 
treated in Ref. p7 while Ref. @ gives a detailed analysis 
of the collective emission of photons through spontaneous 
Raman scattering following excitation by free-space light. 
Section lll Bl orovides an extension of the theoretical treat- 
ment of Ref. |6| to better account for our experimental 
conditions. The emphasis here is the modeling of the de- 
coherence process due to external magnetic fields, and in 
particular for experiments using magneto-optical traps. 
To model this decoherence, the essential elements to be 
introduced in the previous theory of Ref. are the Zee- 
man structure of all levels and an explicit treatment of 
the reading process. On the other hand, the theory in 
Sec. Ill Bl is a simplification of the treatment of Ref. [f| 
concerning the spatial mode of the photons. We consider 
only the forward, collectively enhanced emission. The 
reading process is also treated in a simplified, perturba- 
tive way, while the experiments are done with stronger 
read pulses on resonance. This later difference between 
theory and experiment will result in some noticeable dis- 
crepancy in Sec. IIII B 21 where we discuss measurements 
of the two-photon wavepacket of the pair-generation pro- 
cess. In general, however, the comparison between the- 
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ory and experiment performed in Sec. IIII Bl results in 
very good agreement, which indicates that the theory 
in Sec. Ill Bl takes into account the essential physical ele- 
ments behind the decoherence process. 

A. Photon pair generation 

The building block of the DLCZ protocol is an ensem- 
ble of N identical atoms with lambda-type energy level 
configuration as shown in Fig. ^ which we briefly discuss 
here in an ideal setting. In the experiments discussed 
in this article, the lower states \g) and \s) are hyperfme 
sublevels of the electronic ground state of Cesium atoms. 
First, all atoms are prepared in the state \g). By sending 
in a weak, off-resonant laser pulse, one atom of the en- 
semble might be transfered from \g) to \s), thus emitting 
a photon (field I) at a frequency or polarization differ- 
ent from the original exciting field. A key element of 
the protocol is the collective enhancement of this spon- 
taneous Raman scattering in a forward direction, which 
is determined by the spatial mode of the laser pulse and 
the geometry of the excitation region |(J. If the laser 
intensity is low enough so that two excitations are very 
unlikely, the detection of the photon generated in this 
process is a signature that the ensemble was excited to 
a symmetrical collective state 0, |(| , which in the ideal 
case can be explicitly written as 

1 N 

|la> = -/=$>>!••• ^•••bK, (1) 

where the sum goes over all atoms addressed by the laser 
pulse, and |f a ) indicates the state of the atomic ensemble 
with just one excitation. This is the "writing" step of the 
protocol (Fig.nji). 

Since the excitation probability \ is very small, the 
whole state of the system consisting of atoms and 
forward-scattered mode of light is in the following form: 

|0) = |O a )|Oi) + e^Vx|la)|li) + O(x), (2) 

where \ << 1) \ n i) stands for the state of the forward- 
propagating light field I with n photons, (3 is a phase 
set by propagation to and from the ensemble, and |0 a ) = 
0^° 0(x) represents all the other possible excita- 
tion processes, which in the ideal case occur with proba- 
bilities of order \ 2 - The system remains in this state for 
a time on the order of the lifetime of the ground states. 
By sending in a second ( "read" ) pulse resonant with the 
\s) — > 1 6) transition, the state of the atomic ensemble 
can be transferred deterministically (read out) to another 
forward-propagating light field 2 at the \b) — > \g) tran- 
sition (see Fig. [Qa). In this way, it is possible to access 
the quantum state of the atoms. This reading process 
is then closely related to lo w-lig ht-level Electromagneti- 
cally Induced Transparency |25tl26j|. After the read out, 




FIG. 1: Relevant level structure of the atoms in the ensemble 
for (a) writing and (b) reading processes, with \g) the initial 
ground state and |s) the ground state for storing an excita- 
tion. | a) and \b) are excited states. The transition \g) — > |o) 
is initially coupled by a classical laser pulse (write beam) de- 
tuned from resonance, and the forward-scattered Stokes light 
(field I) comes from the transition \a) — > \s), which has dif- 
ferent polarization or frequency to the write light. A classical 
read pulse then couples the transition \s) — > \b), leading to the 
emission of forward-scattered anti-Stokes light (field 2) from 
the transition \b) — > \g). 

the state of the system becomes: 

|0) = |O 1 )|O 2 )+e^Vx|li)|l2>+O(x), (3) 

where 7 is a phase that includes p and the propagation 
phases to and from the ensemble related to the reading 
process. Fields I and 2 exhibit now strong correlations in 
the photon number basis, and can be described as photon 
pairs. These non-classical correlations can be measured 
by photoelectric detection. Since the field 2 maps the 
state of the atoms, the correlations between field I and 
field 2 can then be used to infer correlations between field 
I and the collective atomic excitations in the sample. 

B. Decoherence 

In order to analyze the decoherence process in the gen- 
eration of pairs from an atomic ensemble as described in 
Sec. Ill Al we need to expand the theoretical treatment 
of Rcf. 6] to include other experimentally relevant fea- 
tures. For our experiments in particular, it is essential 
to include the splitting of the Zeeman structure of the 
atomic ground states due to the magnetic field. The 
MOT quadrupole field generates an inhomogeneous dis- 
tribution of splittings throughout the ensemble. As the 
system evolves in time, this results in dephasing between 
different regions of the atomic cloud, and in a respective 
decay of the coherence of the collective state. It is also 
important to include explicitly the reading process in the 
theory. For simplicity, this is done by considering a read 
process similar to the write process, i.e., with small prob- 
ability of excitation and detuned from the excited state. 
Note that in the actual experiment, the read beam is 
stronger than the write beam and is on resonance. This 
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will lead to small discrepancies when comparing the ex- 
perimental results to the theory, that will be discussed in 
section UTTbH 

The inclusion of Zccman structure in the theory al- 
lows a detailed discussion of the effect of light polariza- 
tion in the experiment. This is important to evaluate 
different excitation and detection schemes. It also gives 
a better description of the initial state, and of its role 
on the subsequent coherent pair generation. Together, 
the analysis of different polarization schemes and of dif- 
ferent initial states led to specific proposals of ways to 
improve the whole process. These features of the theory 
are not specifically related to the MOT magnetic field, 
and should apply to pair generation in other systems, 
like vapor cells or dipole traps. 

Our treatment starts by considering a sample of N 
four-level atoms, such as in Fig. [5] The four levels rep- 
resent manifolds of Zeeman sublevels and are indicated 
by their respective F quantum numbers. A specific state 
of the Fj manifold of the i-th atom is represented by its 
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Two pumping fields act on the system, namely a write 
field E ga and a read field £ s &, where 



£ga{r,t) = u w (f,t)e 



i(k w z 



Pr ) 



(4a) 



(4b) 



which couple the transitions F g — > F a and F s — > Ff,, 
respectively. The functions u w and u r give the slowly- 
varying envelopes of the write and read pulses, respec- 
tively, and e Pm and e Pr are their polarization vectors. As 
a result of their action, two Raman fields are sponta- 
neously generated in the sample: 



£ sa (r,t) oc 



£ g b(r,t) oc 



pi 



E 

P2 



dkxCbz e^ fel r w *i e v . 
1 kipi p 1 



z K2P2 P2 ' 



(5a) 



(5b) 



where u>r = \ki\c and pi is a label for the field polariza- 
tion, ar and br are the annihilation operators for 

feipi k 2 p2 r 

the Raman fields 1 and 2, respectively, which couple the 
transitions F s — > F a and F g — > F^. The state of field 1 
with just one photon excited in mode k\p\ will be desig- 



nated by |l£ 
2. 



kipx 



A similar notation will be used for field 



The Hamiltonian for the system of N atoms can be 
written as 



H(t) =H + V(t) 



(6) 
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FIG. 2: Energy level scheme considered for the atomic ensem- 
bles 



, where rrij is the azimuthal quantum number. where 



A' 



H = ^2 I ^2 (-hvs + ^Bg s m s B Zi ) \m s )i(m s \ 

i=l I m B = -F s 

+ E VB9gmgB z% \m g )i(m g \ 
m g =-F g 

F a Ft, 

+ huj a \m a )i(m a \ + ^b\mb)i(m},\ 



m a = -F a 



m b = — F b 



(7) 



is the free-atom Hamiltonian, and 



JV 



V(t) = J2\ E E {-dm a m g ■ £ga) \m a ) l (m g \ 

i—1 \ m a — — F a m g — — F g 
F s F a 

+ E E (-dm 3 m a ■ £L) \ m s)t(m a \ 

m s — — F s m a — — F a 
F b F B 

+ ^2 E \-dm h m a ■ £sbj \mb)i{m s \ 



JTi b = — Fb Tn$ — — F s 

E E (- 

m g =—Fg m,b = — Fb 



•2gb) KKM \ (8) 



gives the time-dependent interaction Hamiltonian. djk 
is the dipole moment for the j — ► k transition, [ib the 
Bohr magneton, gj the hyperfine Lande factor for level 
Fj, and B Zi is the magnetic field in the position of the 
i-th atom. The magnetic field direction is taken as the 
quantization z axis. We neglect the Zeeman splitting of 
the excited states since we want to investigate a situation 
where it is always smaller than the excited-states natural 
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linewidths. The factors — djk ■ £kj can also be written as Hamiltonian Vi(t) in the interaction picture as 



-d ■ £ — K V {?■ f)p i { k '» x *- u '» t ) (Qa) 

u "m a m g ^ ga — xv m a m g \> v jc , V / 



Um s m a O sa — ^ / w l A ra 1 iii a »j: ipi e 1 ) 

Pi 

(9b) 



-i(k 2 - r—wz t) 



K2P2 



(9d) 



where if ma m 3 , K^ la , K mbms , and if™ 2 ^ are coupling 
constants for the corresponding transition. 

The temporal evolution of the coupled system consist- 
ing of ensemble + Raman fields is described by the evolu- 
tion of its density matrix p(t). In the interaction picture, 
the corresponding operator pi(t) is given by 



p I {t) = u I {t)mv\®-. 



(10) 



where Ui(t) is the temporal evolution operator, and the 
initial state p(0) can be written as 

p(0) =M(0)®^F a (0)®pi(Q)®pa(0)®---®Pjv(0), (11) 

with £fi(0) the initial state of field 1, pf 2 (0) the initial 
state of field 2, and Pi(Q) the initial state of the i-th atom. 
For most of what follows, we will be interested in the 
case where the fields 1 and 2 are initially vacuum states, 
Afi(0) = \vaCF 1 )(vaCF 1 \ and Pf 2 (0) = \vacF 2 )(vaCF 2 \, 
and all atoms are initially in the same incoherent dis- 
tribution over the Zeeman sublevels of the F g state: 



,(0)- £ D 



m g | ' ' L g / i \" L g \ ; 



(12) 



m„ = ~F„ 



with D m giving the probability of finding an atom in the 



m g state at t = 0. In section lTVl however . we will consider 
the case where all the atoms are optically pumped in one 
of the Zeeman sublevel (ttt-f = 0). 

The operator U (t) can be written as a Dyson series in 
the form 



N 



N 



(13) 



i=l i=l 

where 

U?\t) = J\t' j\t"V l {t')V l (t"), (14) 

and so on. The single- atom interaction operator Vi(t) is 
defined from the expression for the general interaction 



/v 



V r (t) = e iHot/hy^y-iH t/n = ^ p.^ _ ^ 



1. Probability for joint detections 

We want to calculate in the lowest order of perturba- 
tion the probability of detecting a single photon in field 
1 followed by another photon in field 2. The first step is 
then to calculate the restriction of the coupled state pit") 
to the space of states of fields 1 and 2: 



p Fl F 2 (t) =Tl A [p{t)\ 



(16) 



The symbol Tr^ indicates a partial trace over all atomic 
states. The probability for detecting two photons, one in 
mode k\p\ and the other in mode k 2 p 2l up to time t is 
then given by 



= (h 1P1 \(h 2P2 \^A W)]\k 2P2 )\h 1P1 )- 

(17) 



Since all atoms are initially in the ground state F g , the 
lowest order term of series (fT3)) that results in a single 
photon in field 1 and another photon in field 2 is the 
fifth term, which accounts for the four transitions carried 
successively by the write field, photon 1, read field, and 
photon 2, respectively. Substituting Eqs. (|ll)fl and 
into Eq. p7|) and keeping only the lowest order term, we 
arrive then at 



p{ 2 (t, k 1 p 1 ,k 2 p 2 ) = 

N 



(18) 



Note that acts only over the fc-th atom. Thus, the 
trace Tr^ on each term of the double sum can be written 
as a trace Tr^ over the states of the atoms at which 
the i/^ operator is acting, since all other atoms remain 
in their initial state. Two different cases are present in 
Eq. 11811 . If i ^ j, the two operators act over two different 
atoms and the initial state p(0) simplifies to PFi(O) <S> 
Pf 2 (0) ® MO) ® Aj(0). If i = j, then p(0) -> p Fl (0) <g> 
Pf 2 (0) 8> Pi(0). With these observations in mind, we see 
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that Eq. 118(1 can then be written as 
p t 1 2(t,k 1 p ll k 2 p 2 ) = 

N 

E (h lPl \(h 2P2 \ T ^ [# 4) (*)M0)] \vac F2 )\vac Fl ) 

C4)t/ 



N 



^(l gipi I (1^ |Tr, W| 4) (<)/5 Fl (0) <g> M (0) ® ft(0) 



x U^Ht)] \l t )|L- ). 
Substituting Eq. ljT2"|l . we have 



(19) 



N 



E E ^mjAifmj.mj) 



JY 



i=l m g = -F g 



EE E A^K.rn,)! 1 



i— 1 m'— — F Q m Q ——F a 



N 

E 



E D mg Ai(mg,m g ] 



m„= — F„ 



(20) 



where 



A % {m' g ,m, g ) = 

(lfe 1P J<lfe P JKI# ) (*)l m s)^ ac ^)l Mc ^)' ( 21 ) 

Note that the first term on the right side of Eq. i|20|) 
scales as TV 2 , while the two remaining terms scale with 
N only. Since we are interested in the limit of large N, 
we can then approximate 



N 



E E D mg A i( m g> m g) 



i=l m„=—F„ 



(22) 



Thus, for large N, only transitions that start and end in 
the same state contribute to the pair generation. This 
result can be understood as a constructive interference 
between all pathways that connect the ensemble back to 
its initial state, after which it is not possible to distin- 
guish which atom made the transition |27j . Pathways 
connecting different initial and final states leave a trace 
in the ensemble, which in principle can give information 
on which specific atom made the transition. In this last 
case, the number of possible pathways generating the pair 
of photons is then linearly proportional to the number of 
atoms N. Eq (|22|l expresses the collective enhancement 
that is essential to the scheme of ref. 01 • 



Finally, substituting the specific expressions for (t) 
and Vi(t), we find that Ai(m g ,m g ) can be written as 



A i (m g ,m g )= 



d{m g ,ms) i(fc) 



x / at e fc 2 r ^ zg > 
Jo 

x [ dt"u r (f f i ,t")e i(Ar - au)t " 
Jo 

x / rf i"' e i(%r A «+»")'"' 
Jo 

x /' ^(f/)^-^, (23) 
Jo 

where A„, = u) a — u> w , A r = cjj + lu s — lj Ti Acj^ = 
u>£ ~lo w -uj Si Auj£ 2 = uj l t 2 - 0J r + uj Si and 

F b F a 

d(m g ,m s ) = Krngrn b Km b m B Km^m a Km a m g 
m h — — F h m a =—F a 

(24) 

gives the strength of an specific excitation pathway in 
which the atom starts at m gi then goes to m s , and 
ends at m g again. The Zeeman splittings are written 
in terms of the parameters ai g = ^BggTn g B Zi /h and 
ai S = fiBg s m s B z Jh. 



2. Forward emission 

In order to simplify the following analysis while keep- 
ing the essential trends of the temporal dynamics, we 
will focus now on the treatment of the forward, reso- 
nant emission from the atomic ensemble. In the forward 
direction, the light emitted by the sample satisfies the 
phase-matching condition 



(25) 



The resonant conditions for the Raman fields are Awr — 
and Aw£ 2 =0. A discussion about deviations from 
these conditions can be found at Ref. 

Under these assumptions, and with the slow envelope 
functions written as 



u r (fi,t) = q r {n)fr{t) , 
u w (n,t) = q w (fi)f w (t) . 



(26a) 
(26b) 



Equation (|2*5|) becomes 



Ai{m g ,rn g ) = q r [r t )q w [r t ) 2^ F(t,Zi), 



m s — — F s 



(27) 
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with 

F(t,Zi)= fdt'e i( -- Ar+ai ^ t ' f dt"/ r (t")e i(Ar ~°** ) *" 
Jo Jo 

x / < ft'V(-A«+a«.)t'" f dt' 1 ' f w {t lv )e i ^- a ^ t ' u . 
Jo Jo 

(28) 

Note that the F function depends on the parameters for 
a specific atom only through Zi that specifies its position 
along the quantization axis. In this way, after a certain 
time, atoms in different parts of the ensemble contribute 
to the probability amplitude of the process with different 
phases. 

If we consider a uniform distribution of atoms through- 
out the beam path, and neglecting the z dependence on 
the q functions, the sum over all atoms may be trans- 
formed in the following integral 



N 



N 



i=l 



dx dy dz q r (x, y)q w (x, y) 



dxdy q r (x,y)g (x,y) N } ^ 

J\ 1 

f L/2 dz 

(q r (x,y)q w (x,y)) N / — 



L/2 



(29) 



where V = AL gives the volume of the excitation region, 
A its transverse area, and L its length. 

Substituting Eqs. JUl and JUJ) in Eq. we finally 
obtain 



pW) = c 



where 



F " Fg f L/2 dz 

V V D m d(m g ,m s ) / — F(t,z) 
^_ Fam X-F, J-L/2L 



C = N 2 \(q r (x,y)q w (x,yW 



(30) 



(31) 



is a constant. After the read pulse has left the sample 
(i.e., when t — > 00), Expression i|3U|l is then proportional 
to the total probability of detecting the pair of photons 
in one trial. Details on how to compare this expression to 
the experimental results will be discussed in Sec. IIII B II 
In the experimentally important case of square pulses, 
it is straightforward to obtain analytical expressions for 
both F(t, z) and pi 2 (t) in the limit of large A w and A r . 



3. Probability density 

Equation (|3l]|) gives the total probability of detecting 
one photon in field 2 after detecting a photon in field 1. 
Now we want to obtain the probability of finding photon 
2 between times t 2 and t 2 + At 2 and photon 1 between 
times t\ and t\ + Ati, for small At 2 and Ati. 



The first step in this calculation is to note that Eq. i|30|) 
can be written as, 



vim = w)v 



(32) 



The function <j)(t) gives then a probability amplitude for 
the process where the two photons are found up to time i. 
It consists of an integral over all possible pairs of detec- 
tion times (^2 , ti ) , representing different excitation path- 
ways, and can in principle also be written as 



4>(t) = dh dt 2 P(t 2> ti) 



(33) 



where we considered explicitly t 2 > t\. P(t2,ti) repre- 
sents then a density of probability amplitude. 

The probability amplitude for finding photon 2 be- 
tween times t 2 and t 2 + At 2 , and photon 1 between times 
t\ and t\ + Ati, can be obtained then by restriction over 
the temporal integral in Eq. (|30fl . Since all the tempo- 
ral dynamics in Eq. is m the function F(t,z), we 
need to calculate first the restriction of F(t, z) for these 
specific processes. In order to do so, note that, in the 
fourth order integral of F(t, z), the emission of photon 2 
is described by the last integral (over t'), while photon I 
emission is described by the third integral (over t'"). The 
restriction of F(t, z) for the emission of photon 2 between 
times t 2 and t 2 + A^2, and photon 1 between times ti and 
h + Ati, is then given by |28| 

G(t 2 ,At 2 ,t 1 ,At 1 )= / dt > e i(-A r+aig )t> 
Jt 2 

x f t dt"f r (t")e^- a ^ t " f tl+Atl dt "' e i(-^+^s)t"' 
Jo Jti 

x / dt' u f w {t lv y { ^"- a ^ )tlv . (34) 
Jo 

Equation l|34l) can be directly evaluated for the case of 
square pulses and large dctunings, such that A r , A w >> 
At^~ ,At± . If the time intervals are also small when 
compared to the timescale of oscillations determined by 
the Zeeman shifts (i.e., At 2 ,Ati << a~Z 1 ,aJ 1 ), then 



Eq. I|34fl can be written as 



G(* 2 , At 2 ,h,Ah) = g(t 2 , h)AhAt 2 , (35) 



with 



In this case, F{t, z) can be derived by 
F(t,z)-- 



t ,t 
dti / dt 2 g(t 2 ,ti) 
Jti 



(37) 



An important remark is that, since any pulse envelope 
can be approximated by a sum of square pulses of differ- 
ent intensities and small duration, Eq. I|36|l is indeed valid 
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for arbitrary pulse shapes, as long as the envelope tem- 
poral variation occurs in a much longer timescale than 
Aii or At 2 . 

The connection between g(t2,t±) and the density of 
probability amplitude P{t2, t\) is then made through the 
relation 

_ Fg Fg 
Pfa,t 1 )=VC J2 J2 D mg d(m g ,m s ) 

trig — — F g m s — — Fs 

r L/2 dz 

x / -r9(t2,ti). (38) 

J-L/2 L 

Finally, the probability density for detecting one photon 
from field 1 at time t\ and another from field 2 at ti is 
associated to 

V(t 2 ,h) = |P(i 2 ,t!)| 2 . (39) 

This is the quantity to be compared with the experimen- 
tal results of Sec. IIII B 21 for the two-photon wavepacket 
of the photon pair. 

III. EXPERIMENTS 

Up to now, the experimental implementation of the 
DLCZ protocol in MOTs have be en p lagued by extremely 
short coherence times ^2 13 Il6l Il7j . As discussed 
above, this short coherence time is a result of the action 
of the MOT quadrupole field over the Zeeman structure 
of the hypcrfinc ground states. In the following, we are 
going to describe a series of experiments that allowed us 
to obtain photon pairs from the trapped atomic cloud 
in a situation of very small magnetic field. In this way, 
we were able to measure coherence times of more than 
10 fxs (more than two orders of magnitude longer than 
the duration of the excitation pulses), and two-photon 
wavepackets for the photon pairs that do not exhibit dis- 
tortion by decoherence even when write and read pulses 
cease overlapping in time ^3] • 

The crucial point is to turn off the MOT magnetic 
field and determine the experimental conditions with a 
best tradeoff between high repetition rate and high op- 
tical density. Note that the atoms fly away from the 
trap and the density starts to decrease when the mag- 
netic field is turned off. Hence, the MOT field has to 
be turned off as fast as possible, to decrease the tran- 
sient time and maximize the region with low magnetic 
field and high density. A fast turning off of the magnetic 
field in our metallic vacuum chambers, however, is not 
straightforward and requires specific techniques, as will 
be discussed in Sec. IIII Al 

Inside each MOT-off period, it is possible to conduct 
many trials of the photon pair experiments. These are 
photon counting measurements that require many events 
in order to acquire good statistics. Hence, we would like 
to have as many MOT-off periods as possible to accumu- 
late a large number of trials. However, the MOT needs 



some time to recover its original density after each off 
period, and this time limits how often it can be turned 
off while still keeping a high enough atomic density. 

During the process of turning off the magnetic field and 
determining the proper conditions for the photon count- 
ing experiments, it was essential to be able to perform 
simpler experiments giving direct access to the ground 
state broadening by the magnetic field. We chose then 
to setup a copropagating stimulated Raman spectroscopy 
apparatus to help us in this process. The results for the 
Raman spectroscopy measurements and the investigation 
to determine the best experimental conditions for the 
photon pair generation are described also in Sec. IIII Al 

The nonclassical correlation experiments are discussed 
in Sec. IIII 51 There we show that the coherence time in- 
creases by more than two orders of magnitude once the 
magnetic field is switched off, and describe measurements 
of the shape of the two-photon wavepacket in both situa- 
tions. In this section, we also compare the experimental 
results with the theory of Sec. Ill Bl 



A. Characterization and magnetic field nulling 

As anticipated above, we use copropagating stimulated 
Raman spectroscopy |29j to probe directly the broaden- 
ing of the hyperfine ground states. Our choice for this 
specific technique is based on the fact that it is insensitive 
to Doppler broadening, but very sensitive to any broad- 
ening caused by magnetic fields, exactly like the spon- 
taneous Raman emission process underlying the photon 
pair generation in our experiment. Raman stimulated 
transitions (see Fig. |3Ji) are two-photon transitions con- 
necting one ground-state hyperfine level to the other one, 
in which a single photon is absorbed from one Raman 
beam and another photon is emitted in the other beam 
by stimulated emission through a virtual level, which is 
located 3 GHz below the Cesium D 2 line in our setup. 

The Raman process is resonant if the frequency differ- 
ence of the two Raman beams equals the ground-state hy- 
perfine interval, around 9.192631770 GHz for Cesium. In 
the absence of collisions and transit broadening, this two- 
photon resonance is very sharp, with a linewidth limited 
only by the power and duration of the Raman beams |29| . 
In this way, since the specific value of the hyperfine inter- 
val for transitions between \m g ) and \m s ) states changes 
with the magnetic field, scanning the frequency of one 
Raman beam with respect to the other gives direct in- 
formation on the frequency distribution of possible two- 
photon resonances dislocated by the magnetic field, i.e., 
on the broadening of the ground state. 

Our setup for Raman spectroscopy is shown in Fig.[3Ji. 
The two Raman beams and a probe beam are coupled 
to the same polarization maintaining fiber, which takes 
the beams close to the MOT and provides good mode- 
matching between them. The probe beam is coupled with 
the same polarization as the Raman field connecting the 
F = 3 ground state to the virtual level, the other Raman 
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field is coupled with the orthogonal polarization. The 
lens at the fiber output focus the beam to a diameter of 
150 /im in the MOT region. After the fiber, the beams 
pass through a 50/50 beam splitter cube. The transmit- 
ted parts of the beams are used as a reference to compen- 
sate for power fluctuations. The reflected part is directed 
to the MOT, forming an angle of about m 3° with the 
quadrupole-field z axis. The shaded area around the z 
axis in Fig. |3t indicates the path of one of our trapping 
beams. The absorption of the probe beam by the atoms 
in the MOT is then measured with a second detector, by 
comparing the probe pulse height with MOT on and off. 



tween the various mj? 



= 0. 




absorption indicates that the Raman pulses succeeded in 
transferring some population from one ground state to 
the other. 

In this way, a plot of the medium optical depth for the 
probe pulse as a function of the detuning between the 
two Raman fields gives a direct measure of the ensemble 
distribution of energies in the ground states. Examples 
of such plots with the MOT magnetic field on and off 
are shown in Figs. Hfc, and EJd, respectively. In Fig. 0Jd 
the Raman pulses are delayed 4 ms from the moment the 
magnetic field was turned off, and the nulling of the field 
was performed using additional bias coils located around 
the MOT and looking for a reduced width of the Ra- 
man trace. From Fig. |3t to 0Jd, the width of the signal 
is then reduced by more than two orders of magnitude, 
from 5 MHz to about 20 kHz. The 20 kHz linewidth of 
Fig. Et>, however, also includes about 10 kHz that comes 
from power broadening by the Raman beams. To mea- 
sure this power broadening, we applied an extra DC field 
in the z direction in order to split the central peak be- 



transitions, and then mea- 



sured the width of the magnetic-field-insensitive transi- 



As mentioned above, the 



FIG. 3: (a) Experimental Raman spectroscopy setup. The 
Raman beams and the probe beam are coupled into a po- 
larization maintaining (PM) fiber and sent trough a beam 
spitter cube (BS). The reflected part is focused into the sam- 
ple with an angle of 3 degrees with respect to the quadrupole 
field z axis, while the transmitted part is used as a reference, 
(b) Relevant level structures and laser frequencies for Raman 
spectroscopy. 

Before the Raman pulses reach the MOT, an optical 
pumping cycle moves the whole atomic population to just 
one of the hyperfine ground states. Note that for the 
following experiments, we make no attempt to optically 
pump the atoms onto a specific Zeeman state. Hence, 
the atomic ensemble is unpolarized and all Zeeman sub- 
states are populated. The action of the Raman pulses, of 
about 150 lis duration and 10 /xW power, then transfers 
some population to the initially empty level if their rela- 
tive detuning matches one of the two-photon transitions 
of the sample. The probe pulse has a duration of 5 /is 
and comes 50 /is after the Raman pulses. It is resonant 
with the cycling transition connecting the initially empty 
ground state to the 6P 3 / 2 level [F = 4 — * F' = 5 if the 
empty ground state is F = 4, F = 3 — * F' = 2 for empty 
F = 3 state] . The probe power is about 50 nW, to guar- 
antee a low saturation of the transition. It is then very 
sensitive to any change in the initial population, and its 
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FIG. 4: (a) Raman trace with the quadrupole MOT mag- 
netic field on. The trace represents the absorbtion of the 
probe pulse following the Raman beams, as a function of the 
Raman detuning S. The line width FWHM is around 5 MHz. 
(b) Raman trace 4 ms after the quadrupole field has been 
switched off. The fitted linewidth is 20 kHz, including 10 kHz 
of power broadening due to the Raman beams 

quadrupole field of the MOT should be switched off as 
fast as possible, in order to maintain the high optical 
density needed for the DLCZ-type experiments. How- 
ever, switching off the magnetic field generated by the 
MOT coils is usually retarded for two reasons. First, 
the current in the coils decays exponentially, with a time 
constant proportional to the inductance of the coils. Sec- 
ond, the field decay time is increased by eddy currents 
in the metallic part of our vacuum chamber. Depending 
on the metallic configuration of chamber and coils, the 
transient period can last for tens of ms. In order to ob- 
tain a faster transient, we use a fast-switching electronic 
circuit [33, 03 ■ This circuit allows a quick reversal of 
the current in the quadrupole coils in order to compen- 
sate for the eddy currents, and resulted in a substantial 
reduction of the transient time in our system. 
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A detailed description of the magnetic field transient 
is given in Fig. [SJi, which plots the Raman scan lincwidth 
as a function of the delay from the moment the field was 
switched off. Figure [SJl then shows the timescale over 
which the ground state has its energy-distribution profile 
changed from Fig.0Ji to Fig.QJs. We can see that after a 
few miliseconds, the linewidth asymptotically reaches a 
plateau, given by the residual DC field in the chamber, 
that we estimate in this case to be on the order of 10 mG. 
The dashed line in Fig. indicates the measured power 
broadening. Shorter transients can be obtained with a 
different metallic chamber configuration (like in Ref. l30f) 
or using non-metallic vacuum chambers. 
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a reasonably small value corresponding to an acceptable 
decoherence time, and the higher limit by the restriction 
that the density should not vary too much throughout 
the region. We accepted a variation of about 30% in the 
density. The lincwidth varies by about 30 kHz in the 
same interval. 

A better cancellation of the magnetic field can in prin- 
ciple lead to even smaller linewidths and, consequently, 
longer coherence times. However, improvements along 
this line will eventually be limited by a different problem: 
the diffusion of atoms out of the excitation region. This 
effect of course depends on the temperature of the sample 
and on the diameter of the excitation beams. In order to 
directly measure this diffusion time, we use again Raman 
spectroscopy. In this case, Raman traces are recorded as 
a function of the delay between the Raman pulses and 
probe. The measurement is done when the magnetic 
field is off, such that there is only one narrow peak in 
the Raman trace, like in Fig.[3Ji. In this case, the area of 
the peak profile is proportional to the number of atoms 
in the excitation region. Figure shows a plot of this 
area as a function of delay. We see that the population 
decays with a time constant of 900 [is, as given by an 
exponential fit to the data (solid line). Note that this 
measurement was done with beams that have 150 fiva. di- 
ameter, while in the correlation measurements described 
later we use beams with 60 /im diameter, leading to a 
diffusion time of the order of 360 /is. 



FIG. 5: (a) Evolution of the ground state linewidth and (b) of 
the optical depth of the sample as a function of the delay from 
the time when the current is switched off in the MOT coils. 
The linewidth is measured with Raman spectroscopy. The 
dashed line represents the measured power broadening due to 
the Raman beams. The OD is determined by measuring the 
absorbtion of a probe pulse in the sample. In both graphs, 
the dashed area represents the window used for measuring 
correlations at the single photon level. 

In order to estimate the optimal region for photon 
counting measurements, it is important to independently 
measure the decay of the optical depth after the magnetic 
field is switched off. In our setup this is done in a straight- 
forward way by turning off the Raman beams and using 
a probe pulse close to resonance with the ground state 
that concentrates all the atomic population. The results 
of such measurement are shown in Fig. [SJj, for which the 
population was initially pumped to F = 4 and the probe 
tuned 10 MHz below the F = 4 — » F' — 5 transition. 
The optical depth measurements in Fig.JSJa were obtained 
from the absorption at 10 MHz detuning and assuming 
a Lorentzian lineshape for the atomic transition with a 
natural linewidth corrected for power broadening by the 
probe beam. 

Together, the results in Figs. and [SJd allow us to 
determine an optimal window for the experiments of 
Sec. IIII 51 i.e., between 3 and 5.5 ms (dashed region in 
both figures) . The lower limit of this region is determined 
by the moment when the residual magnetic field reaches 
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FIG. 6: Diffusion of atoms out of the excitation region. The 
solid line is an exponential fit with a time constant of 900 /is. 
The Raman beam diameter is 150 fj,m. 



B. Nonclassical correlations 

In order to characterize the coherence time of the sys- 
tem for various quantum information applications, e.g. 
for the DLCZ protocol or for generation of conditional 
single photons, the measurements must be performed at 
the single-photon level. In particular, one must know 
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how long a single excitation can be stored in the quan- 
tum memory. For this purpose, we perform correlation 
measurements between fields 1 and 2 as a function of the 
time delay At between write and read pulses, thereby 
probing how the nonclassical character of these correla- 
tions ( and hence of the correlations between field 1 and 
the collective atomic excitations) is preserved during the 
storage process. 

In order to investigate the quantum nature of the cor- 
relations, we use the fact that there exists a well-defined 
border between the classical and quantum domains for 
fields 1 and 2 that can be operationally accessed via co- 
incidence detection, as was first demonstrated in the pi- 
oneering work by Clauser j^. In this way, we measure 
the joint detection probability pi2 for detecting a photon 
in both fields 1 and 2 in the same trial, and the proba- 
bilities pi and p2 to register a single detection event in 
field 1 and field 2, respectively. By splitting field i with 
a 50-50 beamsplitter and directing the output to the two 
detectors, the joint probabilities pa are also measured, 
where i = 1 or 2. Fields for which the Glauber- Sudarshan 
phase-space function is well-behaved (i.e., classical fields) 
are constrained by a Cauchy-Schwarz inequality for the 
various probabilities H3 , namely: 

R=^1<1, (40) 

311 322 

where 511 = pn/pj, 322 = P22/P2, 012(f) = P12/O1P2), 
and t denotes the time separation between the detection 
of photons 1 and 2. In our system, gn — 522 = 2 in the 
ideal case. However, in practice, 311 and 322 are mea- 
sured to be smaller than 2, due to various experimental 
imperfections. Hence in our case measuring g 12 > 2 her- 
alds nonclassical correlations, and in the following we will 
use this quantity as another figure of merit to quantify 
the loss of coherence in the quantum memory. 

The experimental setup used to measure nonclassical 
correlations between fields 1 and 2 is shown in Fig. [7\ As 
already mentioned the sample consists in a cold atomic 
ensemble of Cesium atoms in a magneto-optical trap. 
Each trial consists of a period of cooling and trapping, 
and of a period of measurement during which all the 
beams responsible for cooling and trapping the atoms 
are switched off. During the measurement period, the 
atoms are initially prepared in level \g) (F=4) by optical 
pumping with a laser beam resonant with the transition 
6S 1/2 (F = 3)^6P 3/2 (F' = 4). _ 

A laser pulse with 150 ns duration from the write beam 
then illuminates the sample. The write beam is tuned 
near the \g) — > |o) (corresponding to F = 4 — > F' = 
4 of the £>2 line, at 852 nm) and induces spontaneous 
Raman scattering to the initially empty level |s) (F = 3). 
The intensity of the pulse is made sufficiently weak, such 
that the probability of creating more than one excitation 
in the symmetric collective mode is very low. After a 
variable delay At, the stored excitation is converted into 
a photon in field 2, by sending a read pulse tuned to the 
transition |s) — > \b) (corresponding to F = 3 — > F' = 4 



transition of the D\ line, at 894 nm). The write and 
read beams are orthogonally polarized and combined at 
the polarizing beam splitter PBS 1 (see Fig. 0). At PBS 
1, the write and read beams are spatially mode-matched 
with a measured overlap of about 93%. The beams are 
focussed to a waist of about 30/xto in the sample region. 




FIG. 7: Experimental setup. Write and read pulses propagate 
sequentially into a cloud of cold Cs atoms (MOT), generating 
pairs of correlated output photons 1 and 2. The write and 
read pulses have orthogonal polarizations, are combined at 
polarizing beam splitter PBS1, and then focused in the Cs 
MOT with a waist of approximately 30 \im. The output fields 
are split by PBS2, which also serves as a first stage of filtering 
the (write, read) beams from the (1,2) fields. For example, 
field 2 is transmitted by PBS2 to be subsequently registered 
by detector D3 or D4 while the read pulse itself is reflected 
at PBS2. Further filtering is achieved by passing each of the 
outputs from PBS2 through separate frequency filters. SM 
stands for single mode. 

After the MOT, fields 1 and 2 are detected at the 
two different outputs of PBS 2. A challenging aspect 
of the experiment is to separate the classical pulses from 
the weak nonclassical fields, since they are temporally 
and spatially overlapped, and their frequencies are only 
9 GHz apart. This is done in several steps, which are 
explained in detail in Refs. and 0. After 

the filters, fields 1 and 2 are coupled into optical fibers, 
split by 50/50 fiber beam splitters, and detected by four 
single-photon Silicon avalanche photodiodes (APD). Fi- 
nally, the electronic signals of the APDs are sent to a 
data acquisition card, in order to record the detection 
events and analyze the correlations. 



1. Coherence time measurements 

In order to characterize the system's coherence time, 
we measure 312 and R as a function of the delay At be- 
tween write and read pulses. We then compare the the- 
oretical quantity p~\2{At) = ^{^(At) to the measured 
312 (Ai) by way of a single overall scaling parameter £ 
for all At, as the rate of single counts in fields 1 and 2 
(pi and P2) is measured not to depend on At, to within 
20%. In Fig. [SJl we show our results for 312 with the 
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MOT magnetic field on together with the corresponding 
theoretical fitting. This figure was presented already in a 
previous article [l3| and shows a fast decay of the coher- 
ence between fields 1 and 2, taking place in a time scale 
of less than 200 ns. Note, however, that the coherence 
time is actually smaller than 100 ns, since the write pulse 
itself has a duration of 150 ns. The repetition rate of the 
trials in this case is 250 kHz. The rate of coincidence 
events (detection of photon 1 and photon 2 within the 
same trial) is between 2 and 3 counts per second. 

The theoretical joint probability p*2 i s calculated from 
Eq. (|30|) . assuming C = 1. In this way, we need to per- 
form integrals of the F function over the z coordinate. 
This function depends on z only through the parameters 
a g and a s . The atomic ensemble is assumed to be ini- 
tially unpolarized, i.e., with the atoms evenly distributed 
among all Zeeman states of the \g) level. For the ground 
states of Cesium, we have that the hyperfine Lande fac- 
tors g g and g B of levels \g) and |s), respectively, are given 
by fisgg/h = —fiBgs/h = 0.35MHz/G, so that we can 
write 



a g = l-KKrrig 



(!) 



a s = —2irKm s \ — 
- L 



(41a) 
(41b) 



where we considered the magnetic field for the MOT in 
the form B z — bz, with b the field gradient in the center 
of the MOT, and the constant K given by 



K = 



VBg g 



bL 



(42) 



The value of Kmp gives an estimate for the inhomoge- 
neous broadening associated with level \F, mp) due to the 
magnetic-field gradient b. Note that writing a g and a s as 
in Eqs. (|41|l allows us to perform all spatial integrations 
over the dimensionless coordinate s = z/L, and to com- 
bine many of the relevant experimental parameters in a 
single parameter (K). For our experiment, L = 3.6 mm 
and b = 8.7 G/cm, so that K = 1.1 MHz. This K value 
is consistent with the measurement of the ground-state 
broadening shown in Fig. 

The solid curve in Fig. |HJt shows the theoretical fit- 
ting of pi2(Ai) to the experimental data. We considered 
if = 1.1 MHz in the theory, as estimated above for our 
experimental conditions. The only fitting parameter used 
was £, which was found to be £ = 1.05 x 10 s . Note that 
the theoretical quantity pf^ gives the probability for joint 
detection of the two photons, while 512 is a measure of 
this joint probability normalized by the probability of un- 
correlated coincidence detections. Thus the scaling factor 
£ should be given roughly by the inverse of the probability 
for these uncorrelated coincidences. A theoretical estima- 
tion for this value is given by £ th — [p'^(A* — > 00)] _1 , 
i.e., the inverse of the theoretical joint probability af- 
ter the coherence has completely decayed. For the solid 
curve m Fig. Hi, we find £, th = 1.96 x 10 8 . The difference 



between £ and £*" can be attributed to other sources of 
uncorrelated coincidences (such as dark counts in the de- 
tectors, or leakage from the filters) that are not accounted 
by the theory, which leads to £ < £ . It is also impor- 
tant to have in mind that the noise floor is higher when 
the pulses are overlapping, since there is more leakage 
from the filters in this condition. This results in some 
extra discrepancy when comparing theory to experiment 
by means of one single scaling parameter to all regions of 
Fig. Hi. 
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FIG. 8: Measurement of g\2 as a function of the storage 
time, (a) with the quadrupole field on (taken from \u\\ and 
(b) with the quadrupole field off. The observed decay in (b) 
is consistent with the residual magnetic field in the chamber, 
as measured by Raman spectroscopy. 

The <7i2(Af) measurements with magnetic field off are 
presented in Fig. 03. In this case, we use the information 
acquired from the investigation of Sec. IIII Al and turn off 
the field for a duration of 5.5 ms, at 40 Hz repetition rate. 
From the magnetic-ficld-off period, we use for correlation 
measurements only the 2.5 ms window shown in Fig. 
This 2.5 ms window is then divided in 208 trial periods 
of 12 /is, which results in an overall repetition rate of 
8.3 kHz. In the beginning of each trial, the trap light of 
the MOT (tuned in the F = 4 to F' = 5 transition of the 
D2 line) is turned on for 0.6 (is, and its repumper laser 
(tuned from F = 3 to F' = 4) for 1 /is. This procedure 
prepares the system in the proper initial state, with all 
atoms at the F — 4 hyperfine level of the ground state. 
In this case, the rate of coincidence counts drops to about 
0.33 coincidences/s. 
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Figure |HJd shows then an increase of more than two 
orders of magnitude on the coherence time of the sys- 
tem, when the magnetic field is turned off. The coher- 
ence time is now limited mainly by the rate at which 
we can turn off the magnetic field, and also to some ex- 
tent by our ability to magnetically isolate the system. 
Note that in Fig.JSJi, the Raman-trace linewidth indicates 
that the magnetic field in the measurement window is 
still decaying. The solid curve in Fig. [HJa gives the decay 
theoretically expected for a magnetic-field gradient such 
that K = 12 kHz, corresponding to magnetic fields of 
the order or smaller than 30 mG acting on the ensemble. 
This gives a reasonable approximation to the behavior of 
<7i2 under the action of the residual magnetic field, even 
though the spatial dependence of this field can be more 
complicated than a simple linear gradient. The change in 
K from 1.1 MHz to 12 kHz is consistent with the reduc- 
tion of the ground state linewidth between the two cases, 
as measured directly by the Raman spectroscopy setup. 
Finally, for Fig. Oft f = 0.67 x 10 8 and £ th = 2.2 x 10 8 . 

From Fig. HJb, we see that the correlations are still 
highly nonclassical after a storage time of 10 fis. How- 
ever, from the theoretical fitting we can infer that 512 
should became smaller than 2 at about 25 /zs, which gives 
an estimation for our quantum memory time. 

As discussed above, the measurements with gi2 > 2 
give a strong indication of the nonclassical correlations 
observed in our system, based on reasonable assumptions 
for <jii and (722- The most appropriate verification of the 
nonclassical nature of fields 1 and 2, however, is given by 
the measurement of R as defined in Eq. I|40|) . Such mea- 
surements with the magnetic field off are shown in Fig.0 
More specifically, in Fig. we show the measurements 
of 511 and §22 for the same data points of Fig. [SJd. Sub- 
stituting the results of Figs. and Eli in J3DJ, we then 
obtain the values of R shown in Fig. 03, which confirm 
the strong nonclassical correlation present in our system 
for more than 10 /is. 

The R measurement presents considerably larger error 
bars than for gi2- This comes from the large statistical 
uncertainties involved in the determination of 322, which 
requires measurement of the two-photon component of 
field 2 [l4j]. For this reason, we decided to carry out a 
much longer run of the experiment for the longest coher- 
ence time we were able to probe, 10 /is, which resulted 
in the considerably smaller statistical error of this point. 



2. Two-photon wavepackets 

Central to the DLCZ protocol is the ability to write 
and read collective spin excitations into and out of an 
atomic ensemble, with efficient conversion of discrete spin 
excitations to single-photon wavepackets. A critical as- 
pect of such wave packets is that they are emitted into 
well defined spatiotemporal modes to enable quantum 
interference between emissions from separate ensembles 
(e.g., for entanglement based quantum cryptography ,5]). 



The high efficiencies achieved in the work of Ref . Q en- 
abled us to investigate in detail the temporal properties 
of the nonclassical correlations between emitted photon 
pairs |l7|. providing a direct look at various important 
features of the two-photon wavepacket (field 1 + field 2) 
generated by the system. In the following analysis, our 
main quantity of interest is p T {t\, £2), the joint probabil- 
ity for photoelectric detection of photon 1 at time t\ and 
photon 2 at time ti within a time window of duration 
r. The times for this quantity are counted starting from 
the beginning of the write pulse. This quantities is de- 
termined from the record of time-stamped detections on 
all four photodetectors. The detectors have a time reso- 
lution of 2 ns (minimum bin size), but usually we need 
to consider larger bins to acquire enough events for the 
statistics. 




u 1 1 , 1 1 1 , 1 , I_J 

2 4 6 8 10 
At (US) 

FIG. 9: (a) Measurement of g\\ (open squares) and (722 (open 
circles) as a function of the storage time, (b) Measurement 
of the coefficient R as a function of the storage time. The big 
statistical errors are mainly due to statistical uncertainties in 
the measurement of gn and #22- The points at 10/is have 
been measured for a much longer time and exhibit smaller 
statistical error. 

In our earlier experiments |l7j . we focused on two 
cases: (I) nearly simultaneous application of write and 
read pulses with offset At = 50 ns shorter than the du- 
ration of either pulse, and (77) consecutive (non over- 
lapping) application of write and read pulses with At = 
200 ns. Results for vAU . to) are presented in Fig.lTUl In 
case (I), Fig. HOfa, shows that p T (ti,t2) peaks along the 
line t2 — t\ = 5t\2 — 50 ns with a width At 12 — 60 ns, in 
correspondence to the delay 5ti2 and duration Ati2 for 
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read-out associated with the transition |s) — > |6) — > \g) 
given an initial transition \g) — > \a) — > \s) [T^ . In case 
(77) with the read pulse launched 200 ns after the write 
pulse, the excitation is "stored" in the atomic ensemble 
until the readout. The production of correlated photon 
pairs should now be distributed along ti ~ At + 8t\2 
with width ~ Ati 2 . Instead, as shown in Fig. 110b . 
Pr{ti,t 2 ) peaks towards the end of the write pulse (i.e., 
ti > 100 ns), and near the beginning of the read pulse 
(i.e., 200 < t 2 < 300 ns). Early events for field 1 lead 
to fewer correlated events for field 2, as p T (ti,t 2 ) decays 
rapidly beyond the line t% — ti = — 175 ns. The 
marked contrast between p T {ti,t2) for At = 50 and 200 
ns results in a diminished ability for the conditional gen- 
eration of single photons from excitation stored within 
the atomic ensemble [l II and, more generally, for the 
implementation of the DLCZ protocol for increasing At. 
The underlying mechanism is again decoherence within 
the ensemble. 

By contrast, when the magnetic field is turned off, this 
distortion in the two-photon wavepacket is eliminated 
due to the extended coherence time. We now observe 
the shape shown in Fig. 110b . The delay in Fig. 1 101 is 
At = 1 ^s. 

The theoretical results corresponding to these three sit- 
uations are shown in frames (b), (d), and (f) of Fig. 1101 
These are plots of Eq. (|39[) averaged over 4 ns time win- 
dows for both t 2 and t\, the same time window used 
for the experimental data. We also considered pulses of 
trapezoidal shape, with 20 ns rising time, and FWHM of 
150 ns for the write pulse and 120 ns for the read pulse. 
These values correspond to the experimental parameters. 
The only effect of both the time window and pulse rising 
time is to smooth the edges of the distribution. Differ- 
ently from the case of integrated probabilities, it is nec- 
essary here to introduce more details in the description 
of the pulse shapes, since the theoretical description for 
this signal predicts that it is directly related to the pulse 
profiles [see Eq. Jj2jJ|]. 

The main point that calls our attention in these figures 
is the fact that the theory offers a reasonable explanation 
for the data from consecutive pulses (At = 200 ns) with 
magnetic field on, but not for overlapping pulses or At = 
1 jits with magnetic field off. This discrepancy can be 
simply understood, however, if we remember that one 
of the main approximations of our theory is to consider 
low intensities for both write and read pulses. At low 
intensities and zero magnetic field, the theory gives a 
small and constant probability for the photon 2 emission 
after photon 1. From Eq. we see that the magnetic 
field introduces different phases for different groups of 
atoms. These different phases are proportional to the 
time difference between the emission of photons 2 and 
1, and result in an overall decay of the probability of 
emission of the second photon over time. In Figs. 110b 
and II Of . however, we see that the predicted decay time is 
much longer than the one inferred from the experimental 
data. 



On the other hand, for the actual experiment, the high 
intensity of the read pulse should lead to a fast emission 
of photon 2 once the atom is transferred to level F s . This 
is consistent with the short duration of correlation Ati^ 2 
in Figs. 110b and 110b . which can be understood as com- 
ing from the fast depletion of the F s state. However, 
this reasoning cannot explain the shape of Fig.llOb. since 
the strong excitation alone should result in a similar fast 
depletion in the beginning of the read pulse for any de- 
tection time of photon 1 (as seen in Fig. HUfel. The good 
comparison between Figs. 110b and HOt i comes from the 
fact that the decay due to the magnetic field takes place 
before the delayed readout process occurs. The shape 
in Fig. 110b is then a convolution of a uniform excitation 
probability over t\ (like in Fig. HOb l with the excitation- 
probability distribution oflTUtl. 




FIG. 10: Theory and experiment for two-photon wavepackets 
Pr(ti,t2)' (a) Measured two-photon wavepackets for the case 
where write and read pulses are overlaped with a delay of 50 
ns, with the quadrupole magnetic field on. (b) Theoretical 
predictions for the same conditions as in (a), (c) Measured 
two-photon wavepackets for the case of consecutive (non over- 
lapping) write and read puses with a delay of 200 ns, with 
quadrupole field on. (d) Theoretical predictions for the same 
conditions as in (c). (e) Measured two-photon wavepackets 
for nonoverlapping write and read pulses, with quadrupole 
field off. The delay between write and read pulses is 1 /is. (f) 
Theoretical predictions for the same conditions as in (e). The 
vertical scales are given in arbitrary units proportional to the 
joint probability of detecting photons 1 and 2. See text for 
further details. 
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IV. OPTICAL PUMPING 

The theory developed to explain the data in Fig. |S| 
can also be used to devise new ways to improve the 
system. The inclusion of the Zeeman structure in the 
theory, for example, allows the study of different polar- 
ization schemes for both classical excitation and photon 
detection. It also allows the investigation of the role of 
the atomic initial state on the measured correlations. In 
Fig. ^2 we gi ve two examples of possible ways to im- 
prove the system. The solid and dashed lines in the fig- 
ure represent the two experimental conditions of Fig. |H1 
(initially unpolarized samples with K — 1.1 MHz and 
K = 12 kHz), but now with the same scaling factor. 
The dash-dotted curve shows how the K = 12 kHz curve 
changes if the system is initially spin polarized, with all 
atoms in the \F — 4, rap = 0) state. Note that in this 
case the value of p\ t 2 considerably increases, and the sys- 
tem develops a plateau coming from the predominant 
transition \F = 4, mp = 0) -> \F = 3,m F = 0) -> \F = 
A,mp = 0), which is magnetic-field insensitive. Further- 
more, it is possible to devise a polarization scheme of ex- 
citation that allows only this specific transition for any 
At, e.g. as when the write pulse and field- 1 detection are 
<7 + polarized, and the read pulse and field-2 detection are 
<7~ . This is the case for the dotted curve in Fig. ^] 
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FIG. 11: Variation of pi,2 with the delay At between write and 
read pulses for (solid curve) K = 1.1 MHz and an unpolar- 
ized sample, (dashed curve) K = 12 kHz and an unpolarized 
sample, and (dash-dotted curve) K = 12 kHz and an initially 
spin polarized sample with all atoms in \F = 4, tuf — 0). The 
dotted curve corresponds to an initially spin polarized sample 
classically excited by fields with polarizations such that only a 
magnetic insensitive transition is allowed, see text for details. 
The same arbitrary scaling factor was used for all curves. 

The idealized improvements described by the dotted 
and dash-dotted curves of Fig. ITTI however, will probably 
be limited by two effects which are not taken into account 



by the theory. First, in our experimental setup we should 
see a decay with a timescale on the order of 360 fis due to 
the average time the cold atoms take to cross the 60 /im 
beam diameter of the classical write and read pulses. Sec- 
ond, the theory assumes the presence of a magnetic field 
predominantly in the z direction, which defines the quan- 
tization axis. This can be obtained by a pplying an extra 
DC magnetic field along that direction, [34], |35j but any 
residual transverse field should lead to some decay of the 
plateau. In spite of these restrictions, however, we be- 
lieve that such improvements could lead to an increase 
of more than an order of magnitude over the largest ex- 
perimental decoherence time of Fig. [5] It is also clear 
that there is a benefit in the careful preparation of the 
initial state for the magnitude of the measured correla- 
tions. This is an important point that should also be 
taken into account when considering the implementation 
of the DLCZ protocol in vapor cells. 



V. CONCLUSION 

We have presented a detailed study of the decoher- 
ence processes in the generation of photon pairs from 
atomic ensembles, via the DLCZ protocol of ref . @. We 
have identified the main cause of decoherence for cold 
atoms in magneto-optical traps as being the inhomoge- 
neous broadening of the hyperfine ground states due to 
the quadrupole magnetic field used to trap the atoms. A 
detailed theory has been developed to model this effect. 
We also reported a series of measurement to characterize 
and control the decoherence using copropagating stim- 
ulated Raman scattering. These measurement allowed 
us to switch off the quadrupole magnetic field in a con- 
trolled way. With the magnetic field off, we observed 
highly nonclassical correlations between the two emitted 
photons, for a storage time of up to 10 /is, an improve- 
ment of more than two orders of magnitude compared to 
previous results with cold atoms. Furthermore, contrary 
to all related experiments reported up to now, the coher- 
ence time is now two orders of magnitude larger than the 
excitation pulses duration. This is a crucial step in order 
to use atomic ensembles as a quantum memory to store 
conditional single photon states or entanglement between 
two distant ensembles. 
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